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We are concerned with the optimal control problem of the well known nonlocal thermistor problem, i.e., in 
studying the heat transfer in the resistor device whose electrical conductivity is strongly dependent on the 
temperature. Existence of an optimal control is proved. The optimality system consisting of the state system 
coupled with adjoint equations is derived, together with a characterization of the optimal control. Uniqueness 
of solution to the optimality system, and therefore the uniqueness of the optimal control, is established. The 
last part is devoted to numerical simulations. 
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1 Introduction 

Let f2 be a bounded domain in M. N with a sufficiently smooth boundary dfl, and let Qt = 
Q, x (0, T). In this work we are interested to study an optimal control problem to the following 
nonlocal parabolic boundary value problem: 

du (I) 
— = -Pu , on S T = dUx (0, T) , W 
av 



u(0) = uq , in fi, 



where A is the Laplacian with respect to the spacial variables, / is supposed to be a smooth 
function prescribed below, and T a fixed positive real. Here v denotes the outward unit normal 
and = zaV is the normal derivative on d£l. Such problems arise in many applications, for 
instance, in studying the heat transfer in a resistor device whose electrical conductivity / is 
strongly dependent on the temperature u. The equation ([T]) describes the diffusion of the tem- 
perature with the presence of a nonlocal term. Constant A is a dimensionless parameter, which 
can be identified with the square of the applied potential difference at the ends of the conductor. 
Function f3 is the positive thermal transfer coefficient, which can depend only in spatial variables 
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x or time t, but for the sake of generality we take j3 depending in both x and t. The given value 
uq is the initial condition for temperature. Boundary conditions are derived from Newton cooling 
law, sometimes called Robin conditions or third type boundary conditions. In the particular case 
when f3 = 0, we obtain an homogeneous Neumann condition or an adiabatic condition. Other 
boundary conditions appear naturally, but for the sake of simplicity we consider in this paper 
mixed conditions only. Recall that under restrictive conditions, ([1]) is obtained by reducing the 
elliptic-parabolic system of partial differential equations modelling the so-called thermistor: 



Ut = V.(k(u)Vu) + a(u)\V<p\ 2 , 
V{a{u)V(p) = 0, 



(2) 



where u represents the temperature generated by the electric current flowing through a conduc- 
tor, ip the electric potential, and a(u) an d k(u) the e l ectric and ther mal conductivities, respec- 
tively. For more description, we refer to ( Lacev 1995 . Tzanetis 20021 ). A throughout discussion 
about the history of thermistors, an d more detailed accounts of their advanta g es and applica - 
tions to industry, can be fou nd in (iMadenl ll979j, lShi_el_a]J [l99ji iKwokl 1 19951 . ICimattil HblH ). 
Since the paper of Rodriguesl ( 19921 ). which apparently was the first wh o proved the existence 
of we ak solutions to the system (|2|), several results were obtained. In ( Antontsev and Chipotl 
19941 1 existence and regularity of weak solutions to the thermistor problem were established. We 
remember that existence and unique ness of solution to CD under hyp otheses (H1)-(H3) below 
(cf. Sec. [2|) has been establ ished in ([ El Hachimi and Sidi Ammi 120051). For mo re on existence 
and uniqueness we refer to (|Sidi Ammill2010l . IZhou and Liull201ol . ICimattill201lh . 

Optimal control of problems governed by partial differential equations is a fertile field of re- 
search and a source of many challenging mathematical i ssues and interesting applications ( Lionsl 
197ll . lArantes and Munoz Riveral l201oi . iTroltzschl l2010h . Among essential points in the theory 
we mention: (i) existence, regularity, and uniqueness of the optimal control problem; (ii) nec- 
essary optimality conditions, which consist of the equation under consideration and an adjoint 
system. Existence and reg ularity theory of elliptic and parabolic equations was developed since 



gu 

(jLadyzenskaya et al.lll97ll ). Optimal control theory for the system ([2]) received recently an im- 
portant increase of interest. Results for ([I]) are, however, scarcer and underdeveloped. To the 
best of the author's kn owledge, known result s on the optimal control of a thermistor problem 
reduce to the ones of ( Lee and Shilkinl 1200.51 ). where the term source is taken to be the con- 
trol. In (ICimattill2007h the problem of finding the optimal difference of applied potential to the 
thermistor problem ([2]), in the sens e of minimizin g a suitable cost functional involving the tem- 
perature, is studied. Main result of ( Cimatti 20071 ) gives the optimal system in the simplest case 
of a constant electric conductivity. In addition, a theorem of existence of the opti mal solution 
is giy en in the general case of conductivities depending on the temperature. Paper (jSidi Ammi 
20071 ) investigate s a parabolic-elliptic system similar to ([2]), assuming a particular structure of 
the controls. In (jHrvnkiv et al.ll2008l ). authors considered the optimal control of a two dimen- 
sional steady state thermistor problem. An op timal control p roble m of a two dimensional time 
dependent thermistor system is considered in ( Hrvnkiv 20091 ). In ( Sidi Ammi and Torres! 20071 ) 
a similar problem to ([2]) is studied, consisting of nonlinear partial differential equations resulting 
from the traditional modelling of oil enginee ring within the framework o f the mechanics of a 
continuous medium. The main technique of (jSidi Ammi and Torres! 120071 ) is the adjoint state 
and disturbance method to derive the necessary optimality conditions. Recently, the authors in 
( Homberg et al. 2009/10l ) investigated the state-constrained optim al control of th e thermistor 



problem with the restriction to two-dimensional domains, while in (|Cimattill201ll ) some appli- 



cations to the thermistor problem, and to certain problems of filtration of fluids in a porous 
medium in the presence of the so-called Soret-Dufour effect, are given. However, we are not 
aware of any work or study about the optimal control of (HJ. 

It is known t hat large temperature gradients may cause a thermistor to crack. Numerica l 
experiments in dFowler et al.lll99l IZhou and Westbrooklll997l . iNikolopoul os and Zouraris 2008) 
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show that low values of the heat transfer coefficient j3 results in small temperature variations. On 
the other hand, low values of the heat transfer coefficient leads to high operating temperatures 
of a thermistor, which is undesirable from the point of view of applications. This motivates 
the choice of the heat transfer coefficient as the control, and to consider the optimal control 
problem of minimizing the heat transfer coefficient while keeping the operating temperature of 
the thermistor not too high. 



2 Outline of the paper and Hypotheses 

We consider an optimal control problem with the partial differential equations ([I]): 

(i) The control j3 belongs to the set of admissible controls 

U m = {P G L°°(n x (0, T)) , < m < (3 < M} . 

(ii) The goal is to minimize a cost functional J(/3) defined in terms of u{f3) and (3 as 

J(j3) = / udxdt + / P 2 dsdt . 

J Qt J St 

More precisely, we intend to find /3 € Um such that 

J(P) = min J(/3). (3) 

In Section [3l existence and regularity of the optimal control are established through a mini- 
mizing sequence argument. The energy estimates, in an appropriate space, and then the class of 
weak solutions obtained, allow us to study, in Section HI the optimal control problem and to de- 
rive the optimality system. The obtained necessary optimality conditions consist of the original 
state parabolic equation ([1]) coupled with the adjoint equations together with a characterization 
of the optimal control. In general terms, the approach used here is close to the method used in 
( Hrvnkivl 20091 ) for investigation of the time dependent thermistor problem. Since our objective 



functional depends on u, it is differentiated with respect to the control. We calculate the Gateaux 
derivative of J with respect to f3 in the direction / at the minimizer control f3. We also need 
to differentiate u with respect to the control (3. The difference quotient (u(f3 + el) — u((3)) je is 
proved to converge weakly in H 1 ^) to ifj. As a result, the function ip verifies a linear PDE which 
gives the adjoint system, and an explicit form of the optimal control is determined. Section [5] is 
devoted to the uniqueness of the solution to the optimality system, and therefore the uniqueness 
of the optimal control. Finally, in Section [6] we solve the optimality system numerically for a 
constant case of the optimization parameter. 

In the sequel we shall assume the following assumptions: 

(HI) / : R — > R is a positive Lipshitzian continuous function. 

(H2) There exist positive constants c and a such that c < /(£) < c|£| a+1 + c for all £ € BL 
(H3) u G L°°(n). 

We say that u is a weak solution to (H|) if 

/ ^-vdx + f VuVvdx + f /3uvds = - „ \ f f(u)vdx , (4) 
Jn ot Jn J dn (J n /(it) dx) 1 J Q 

for all v € H 1 ^). We use the standard notation for Sobolev spaces. We denote || ■ ||i>(n) = || ■ \\ P 
for each p E [1, oo]. Along the text constants c are generic, and may change at each occurrence. 
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3 Existence of an optimal control 

The proof of existence of an optimal control (Theorem 13. ip is done using proper estimates 
(Lemma 13. 2\) . 

Theorem 3.1: Assume that the assumptions (H1)-(H3) hold. Then, there exists at least an 
optimal solution j3 € L°°(Qt) of ([3]). Function u = u((3) verifies in the sense of distributions, 
with the following regularity: u £ C(0, T, L 2 (fl)), ff € L 2 (0,T, H' 1 ^)), u G L 2 (0,T, H\n)). 

Proof Let ((3 n ) n be a minimizing sequence of J(f3) in Um- In other words, we have 

lim J(/3 n ) = inf J(/3) . 

n-y+oo P&J M 

In order to continue the proof we proceed with the derivation of a priori estimates: 

Lemma 3.2: Let u n = u{(3 n ) be the corresponding solutions to the weak formulation of ([Tj). 
Then ||^ti II^^q t J ff 1 (0)) ll^nlll — c > where c is a constant independent of n. 

Proof Multiplying the corresponding equations of ([I]) by u n and using the fact that u n £ L°°(J1), 
we obtain 



1 — 1 

2dV 



u n\\2 + \ \Vu n \ 2 dx+ / /3 n u 2 ds<c / \f(u n )u n \dx 
Jn Jan Jn 

+ c) \u n \dx 



< 


c 


/( 








Jn 




< 


c 


U n 


1 


< 


c 


Un 


2 


< 


c 


U n 





Using the fact that < m < /3 n , we have 

+ / \Vu n \ 2 dx + m u 2 n ds < c\\u n \\ H Un)- (5) 
Jn Jan 



ld_ 

2di' 



Now denote 



Ml 2 



/ \Vv\ 2 dx + m / v 2 ds. (6) 
Jn Jan 



It is well know n that \\v\\* defines a norm on H 1 ^!) which is equivalent to the || • ||/p(q) norm 
(Zeidler 1988). Then, there exists a constant [i > such that 

II II 2 ^ || ||2 ^ II ||2 (<-,\ 

MlPnllifi(n) S \\ u n\U ^ c \\ u n\\H^{ny \'J 

It follows from ©-([7]) that 

Id., no n2 1 ^ n2 n2 i /U n 

We obtain H^nWIli + /^ll^nll^g th 1 ^)) — c integrating over (0,T). □ 
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We now continue the proof of Theorem 13.11 By Lemma 13.21 we have, for all n, that 



Therefore, f rom (pQ), is bounded in L 2 (0,T,H Using compacity arguments of Lions 



(Lions 1969) and Aubin's lemma, we have that (u n ) is compact in L 2 (Qt). Hence we can extract 
from (u n ) a subsequence, not relabeled, and there exists j3 G Um such that 

u n — y u weakly in L 2 (0, T, H 1 ^)), 
du n du r , 



oi . )f weaJdyhiL 2 (0,r,H- 1 (n)), 

u n — >■ u strongly in L 2 (Qt), 

u n u a.e. in L 2 (Qt), 
(3 n -> (3 weakly in L 2 (dU), 
f3 n ^ weakly star in L°°(dtt). 



(8) 



Our task consists now to prove that u = u(/3) is a weak solution of (JT]) with control /?. From the 
weak formulation of u n we have 



Sit f f X 

-^-vdx + / Vu n Vvdx + / i3 n u n vds = „ — — — , 



f(u n )vdx. 



We first show that for any test function v G i/ 1 (ri) and n — > oo we have 



f3 n u n vds — >■ / /3uvds. 
an Jan 



Indeed, 



j3 n u n vds — I j3uvds 
an Jan 



< 



on 



(/3 n u n v - /3 n uv)ds 



+ 



an 



(f3 n uv — j3uv)ds 



<M [ \u n -u\\v\ds+ [ (fa-P) 
Jan Jan 



)uvds 



< M\\u n - u\\ L 2( dn) \\v\\ L 2( m) \ + 

< M\\u n - u\\ H i( n) \\v\\ L 2( dn) \ + 



on 



{fin — j3)uvds 



on 



fi n — (3)uvds 



(9) 



where we used here the trace inequality ||w||i, 2 (dn) < c||it||#i(n), which gives that u E if 1 (J)) 
implies u G L 2 (dQ,). It is obvious from limits © that the right hand side of the above inequality 
© goes to when n — >■ oo. On the other hand, we have u n — > u a.e. in fi x (0,T). Since / is 
continuous, f(u n ) — > f(u) a.e. in L 2 (Q). It follows that 



f{u n )dx ->■ / f{u)dx, 
n Jn 
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and 

f f(u n )vdx -> f f(u)vdx, Vv G H\n). 
Jn Jn 

We conclude that u = u(/3) is a weak solution of ([T]). Using the fact that J{f3) is weak lower 
semicontinuous with respect to the L 2 norm, it follows that the infimum is achieved at f3. □ 



4 Characterization of the optimal control 

To study the optimal control we derive an optimality system consisting of equation ([1]) coupled 
with an adjoint system. Then, in order to obtain necessary conditions for the optimality system, 
we differentiate the cost functional and the temperature u with respect to the control j3. Here, 
besides (H1)-(H3), we further suppose that 

(H4) / is of class C 1 . 

Theorem 4.1 : Assume hypotheses (Hl)-(H^). Then (3 i— >• u{(3) is differentiable in the sense 
that as e — > 



<^a)-<0)^ raealdy . nHl(fll 



for any {3, 1 £ Um such that (/3 + el) G Um for small e. Moreover, tp verifies 

dt (f n f( u ) dx rJ n numx+ (f n f(u)dxy m ^ 

^ + P^ + lu = Q on dQ. 
ov 

The proof of Theorem 14.11 passes by several steps. 



(10) 



4.1 A priori estimates and convergence 

Denote u = u(/3) and u e = u((3 £ ), where f3 e = /3 + el. Before the derivation of the optimality 
system, we need to establish an H 1 norm estimate of Ue ~ u ■ 

Lemma 4.2: We have 



+ 



< c. 



L 2 (o,T,m<si)) 



Proof Subtracting equation dU from the corresponding equation of u e , we have 

d_ f u £ -u \ _ ( u £ -u \ = A (f(u £ ) - f(u)) A / 1 1 \ 

dt\ e ) \ £ ) e U u f{u e )dxY ^ e J{ > {(j n f(u £ ) dx) 2 (/„ f %u) dx) 2 J ' 

(11) 
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Multiplying the equation (jlip by e , we obtain that 



ld_ 

2di 



u E -u 



V( 



ld_ 

2di 

id_ 

2dt 



2 

+ 

2 

2 

+ 

2 

2 



U F — U. 



u F — u 



£ 

Up — U 



V( 



Up — u. 



dx- V 

2 



Up — U \ Up — u 



v ds 







dx ~ I — 

/an 9^ V e 



Up — U \ Up — u 



£ 



ds 



e 



V(^) 



< 



2 

- /(«) u £ - u 



2 dx + [ p( l^-H) 2 dS + [ IU £ ( I ,1, 

Jan V £ J Jen V e 



(/ n /(«)^) 2 
/(«) 



+ A 



Up — u 



e \(f n f(u £ )dx)* (/„/(«) cfa)V ^ 
Since < m < /3, we get 



ld_ 

2dt 

< 







u £ -u 








£ 


2 ./n 



V 



Up — u 



ld_ 

2di 



Up — u 



+ 
2 JQ 



V 



Up — u 



dx- V 

2 

dx + m 



Up — u\ Up — u 



< 



A / f(u e ) - f(u) u £ -u 



U a f(u)dxY 

/(«) 



+ A 



1 



1 



an V £ 



Up-U 



v ds 



ds + lu £ ( — ) ds 



e \(J Q f(u £ )dxr (f n f(u)dxrj> e 
Using the fact that / is Lipschitzian, it follows from the L°° boundedness of u and u £ that 



2di 



Up -u 



+ 

2 JO. 



V( 



Up-U, 



dx + m 



Up — u 



ds + lu £ 
dn \ £ J Jan V 



Up — u 



ds 



< c 



< c 



< c 





'u £ - 




v £ 


u £ 


— u 




£ 




— u 


£ 



dx H — 



\(f Q f(u)dxf-(j Q f( U£ )dxy 



(kf^e)dxY (/„/(«) da:)' 



/(«), 



Up — u 



(f(u)-f(u £ ))dx) / (/(«) + /(« e ))dx 
2 8 \Jn / \in 



Up — u 



+ c 



Up-U 



Up — U 



< c 



Up-U 



+ c 



e 

Up — u 



£ 



Since L 2 (0) C L 1 ^), then 



2<9i 



u £ -u 



2 

+ 

2 



V 



Up — u 



dx + m 



— ) ds + lu £ i — ) ds 

an \ £ / Jan V e 



< c 



Up — u 
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Using the trace inequality ||ii £ ||L 2 (afi) < cll^ell^fn)! we have 



ld_ 

2dt 



Ue — U 



+ 

2 



V 



Ue—U 



dx + m 



u E -u 



an \ e 



ds 



I an 



ds + c 



u F — u 



< c \\u F 



\L 2 (an) 



+ c 



L 2 (an) 



2 

2 

u F — u 



< c IpeIIh 1 ^) 



u F — u 



+ c 



Up-U 



Thus, 



19 

2dt 



Up — U 



Up — u 



< c \\u e \ 



H!(n) 



Up — u 



+ c 



if 1 (ft) 



Up — u 



On the other hand, by the equivalence of || • ||* and || • ||#i, we have for a positive constant c that 



Up — u 



< 



ff^ft) 



Up — u 



It follows from Young's inequality that 



ld_ 

2dt 



Up-U 



+ c 



Up-U 



ff^ft) 



< c ll n ell,ffi(ft) 



Up-U 



+ c 



/p(n) 



Up-U 



< c 



< 



Up — u 



+ c 



Up — u 



ff^ft) 

2 

H x (n) 



Up — u 



+ c 



Up — u 



+ c. 



Therefore, 



Up-U 



+ c 



Up — U 



< C. 



if 1 (ft) 



We get the intended result of Lemma 14.21 integrating this inequality with respect to time. □ 



Using the energy estimates of Lemma 14.21 we have, up to a subsequence of e — > 0, that there 
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exists ip such that 



£ 

"e 



ip weakly in L°°(0, T, L 2 (Q)), 



" ~ " V 1 weakly in L 2 (0, T, H X (Q)), 
<) (u e -u\ . weaklyinL 2( Q TH -i n))j 

4) weakly in L°°(0,T, L 2 (dn)), 



dt\ej dt 

u £ -u 



e 

•2 



/3 £ -> /3 weakly in I/(dfi) as e — > 0, 
/3 e -> p weakly in L°°(fi) as e -> 0. 



4.2 Proo/ o/ Theorem g77] 

We are now ready to derive system f|10|) . We have 



with 



i(^)- + / n n^)^ + / s /(^)- + i i »« 

= j+jj 



(J n f(u £ )dx) 2 Jn 



and 



//:= — ( 2 " 2) / f( u ) v dx. 

(f n f(u £ )dx) (J n f(u)dx) J Jn 



s 



We can write II as follows: 

e U n f(u)dxf(^f(u £ )dxf J n nU)VdX 
_ A f (f(u)-f(u e) ) J n (f(u) + f(u £ ))dx r 

- L e dxx u n nu £ )dxnu{u)dxy]J (u)vdx 

One can show, using weak convergence (|12|) . that 



-2A Lf(u)vdx f ., x 
Unf( u ) dx ) Jn 



In the same manner we have 

1 x 



(/„/(«) <fa) 2 J n 



/ f'(u)ipvdx as e — > 0. 
in 
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Again, from the weak convergence (|12p . we conclude that, as e — > 0, (j 1 3 j) converges to 



^ vdx + I Vip'Vv + I fitpvds + / luvds 



Jn 9t JQ JQQ JQQ 

— 2 A j Q f(u)vdx 
(J n f(u)dxf J v :' ' <.! {l f(<i)<l:rf Jo 

for every v E H 1 ^). In other words, 



f'(u)ipdx + —j — / f'(u)ipvdx 



f dip f f 

/ -p—vdx + / VtpVvdx + / (flip + lu)vds 
Jn &t Jqci 



2\J n f(u)vdx f fl{u)i)dx+ A . r((7j; , ;v/r (J J) 



We can rewrite (|14|) as follows: 



dip I f dip 

-—vdx+ / — Atpvdx + / (— — \- f3ip + lu)vds 
n ot J n J dn dv 



2A f n f(u)vdx f ,, , A , 

' f{u)ipdx + — / f{u)ipvdx. 



{J Q f{u)dxf J n J (J n f(u)dxyj n 
We conclude that ?/> satisfies the system 

a- A, "= (/„/(„) /( " )+ (/„/(«) m " 
^ + /m + = o on an. 

This completes the proof of Theorem 14.11 



4.3 Derivation of the adjoint system 

In order to derive the optimality system and to characterize the optimal control, we introduce 
an adjoint function ip, defined in Qt and enough smooth, and the adjoint operator associated 
with ip. Multiplying the first equation of (fl"0|) by tp and integrating in space and time, we have 

dip f 
— — ■ (fdxdt + / —Aip ■ (pdxdt 

_-2\f n f'(u)ipdx f J QT \f(u)iP<pdxdt 



/ f(u) ( pdxdt+ JQT ; ;V in». (15) 
JQt C o fW dxY 



(J a f(u)dxf J Qt j " (J Q f(u)dxf 
Integrating by parts (|15p with respect to time, and imposing the boundary and initial conditions 



^ + (3ip = on dn x (0,T), (p(T) = , ip(Q) = , 
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we obtain 



ip(T)ip(T)dx - I ip(0)ip(0)dx + / -^f-.4>dxdt+ I -Aip.i/jdxdt 



2A L r /( n )v 9dx * r , , , ^jn T f'{u)ip^dxdt 
1 f (u)tpdx + 



(/„/(«) cfc) 3 (/„/(«) d*) 2 • 

Thus, the function ip satisfies the adjoint system given by 

dy? a -^.f n f(u)ipdx Xf'(u)ip 

-m- Av= UJWW /(u)+ (j n /(«)^ +lmfc 

(16) 

az/ 

y>(0)=0, y>(T) = 0, 

where the 1 appears from differentiation of the integrand of J(/3) with respect to the state u. 

Theorem 4.3 (Existence of solution to the adjoint system): Given an optimal control j3 € Um 
and the corresponding state u, there exists a solution ip to the adjoint system (|16p . 



Proof Follows by the arguments in (jSidi Ammi and Torres! 120071 ). □ 



4.4 Derivation of the optimality system 

Gathering equation ([I]) and the adjoint system (116j) . we obtain the following optimality system: 

ut — i\u — 



(/„/(«) Cte)2' 

dy? a -^fnf{u)ydx \f'{u)<p> 

-m~^ = (/„/(„>*). /( " )+ (/ n / (t ,)^ +lm ^ 

^ + /?« = 0on^x (0,T), (17) 
az/ 

^ + ^ = 0on3Ox (0,T), 

p(o) = o, <p(X) = o, 

u(0) = u . 

Remark 1 : The existence of solution to the optimality system (|17p follows from the existence 
of solution to the state system ([1]) and the adjoint system (|16p . gathered with the existence of 
optimal control. 

We characterize the optimal control with the help of the arguments of ( Hrynkiv 2009). 
Lemma 4.4: The optimal control (3 is explicitly given by 

P(x,t) = min (max (-^-,m) ,m) . (18) 



V V 2 

Proof Because the minimum of the cost functional J is achieved at /3, using (|10p . the convergence 
results (fT2j) . and the second equation of the system (fTTp . we have, for a variation / £ Um with 
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ft + el € Um and e > sufficiently small, that 
J(P + el) - J(/3) 



< lim 




(u{P + el) - u(P)) dxdt + / ((/3 + el) 2 - /3 2 ) ds dt 

u(f3 + el)-u(8) , , „ f a7 , , 
y -^dxdt + 2 / ^Zdsdt 



xfjdxdt + 2 I /3ldsdt= / (ip,<p)(l,0)dxdt + 2 pidsdt 



< / (2/3/ + /ui^) (is dt 

< 1(2/3 + ucp) ds dt. 
Jdnx(0,T) 



Using the arguments and techniques in ( Hrvnkiv 2009) involving choices of the variation function 



I, we have three cases to distinguish, (i) Take the variation £ to have support on the set {x G 9$7 : 
m < /3(x,t) < M}. The variation £(x,t) can be of any sign, therefore we obtain 2/3 + mp = 0, 
whence /3 = (ii) On the set {(x,t) E S T = dQ, x (0, T) : /3(x,t) = M}, the variation must 

satisfy £(x,t) < and therefore we get 2/3 + u<p < 0, implying M = (3(x,t) < —^f. (hi) On 
the set {(x,t) € dfl x (0, T) : (3(x,t) = m}, the variation must satisfy £(x,t) > 0. This implies 
2/3 + mp > and hence m = (3(x) > — Combining cases (i), (ii), and (hi) gives 




■?f<M, 



p = < M, a ~m>M. 



2 

< m. 



This can be written compactly as f)18[) . □ 



4.5 Particular case: a constant heat transfer coefficient 

Let us consider now the case when the heat transfer coefficient /3 is a constant, i.e., when /3 is 
independent of x and t, and 



J 



/ udx + p 2 . (19) 
in 



We need to adjust the parameter /3 € E/fcf in such way that the new form of the functional (|19[) is 
minimized. Then, all the theory of existence of optimal control and derivation of the optimality 
system, that one needs to put into the proofs of the previous sections carries over to this case 
and are simpler. As for the characterization of optimal control, we have: 

Lemma 4.5: The optimal parameter characterization related to (|19p is 

/3 = min ^max mtj, — — J mpds^J , m \ . (20) 



Proof For the characterization of the optimal control we take into account the new expression 
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of the cost functional (fTUj) : 
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e-S-0 e 



1 



\im-\ I u((3 + el)dx + (p + el) 2 - / u(f3)dx - /3 2 



e e 

V>cfo + 2/3/ 
(^,^)(l,0)dx + 2/3/. 

Multiplying the optimality system (|17p with the test function integrating by parts, and using 
([TO]), we find that 

0<llm W + d)-JW = ,/■ ^ + VL 

o° e Jan 

Therefore, 



Z uipds + 2pj > 0. 

Repeating all the steps as those yielding to (fTHj) . we obtain that the optimal parameter /3 is 
characterized by (l20j) . □ 



5 Uniqueness of the optimal control 

The uniqueness of the optimal control is mainly based on the L°° boundedness of u and ip. These 
are quite realistic as sumptions since physical quantities are always bounded. It has been shown 
m dSidi Ammill2010h that u £ L°°(Q,). It remains to establish that (p is also essentially bounded. 



Lemma 5.1: Under hypotheses (H1)-(H4) one has tp £ L°°(f2). 

Proof Multiplying the second equation of (fTTj) . governed by tp, by — ip k+1 for some enough big 
integer k > 2, we have by the L°° estimate of u and Young's inequality that 

^^IMl£+2 + / Vy>v(M* +1 )dz < cM\v m \W\\ L tl m + 4<p\\ k L i? Hn) + Hl£ft (n) 



< C \\uD\\ k+2 + c\\tn\\ k+1 



Then, 

1 d 



k + 2dt 



\<P\\ttl + ( k + 1 ^J Q ^ 2 \ V V?dx < c\\p\\ k L t 2 2{n) + c\\<p\\ k +l 2(n) . 
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Taking into account that the second term of the left hand side is positive, we have 



1 ^|| \\k-\-2 ^ II nfc+2 || n/i,'+I 



k + 2dt llrllk + 2 ^^»L"+ 2 (Q) ^ MIVMifc+a^)- 

Setting yk = |MlL fc + 2 (f2); it follows that 

„ fc+l °Vk , k+2 , fc+1 

y k -QT-< c Vk + c Vk ■ 

In other words, 

dy k . 

-K- < cy k + c. 
at 

By the Gronwall Lemma we have y k = \\<f\\L k + 2 (n) — c i where c are constants independent of k. 
Letting k — > oo, we have ||v||x,«>(-n) < c - D 

Theorem 5.2: If the hypotheses (H1)-(H4) hold, then the solution of the optimality system 
(|17p is unique and, therefore, the optimal control (3 is unique. 

Proof Let u±, tp± and U2, tp 2 be two solutions to the optimality system (fT7|) and f3±, j3 2 be two 
optimal controls. Denote w = u\ — U2 and <p = (f2 — tpi. Upon subtracting and estimating the 
difference between the equations governed by u\ and U2, we have 

9 a A/(m) A/(u 2 ) , . , . 

*™" A " = (/„/(«.)*)' " 0n7W^P =9( " l '" 2) "'' <21) 

where 

/ \ f A/(m) A/(« 2 ) A , 

VCJn/wW J n f( u v dx J 

By using hypotheses (H1)-(H4) and the L°° estimate of U{, i = 1,2, we have fl^ui, u 2 ) € L°°(0). 
Multiplying (j2Tj) by w yields 



ld\\w 



1 2 



2 Si 
Since m < /3 2 , we have 



+ / |Vw| 2 (ix+ / /3 2 |u>| 2 ds + / (/3 2 - (5{)uiwds < c\\w\\l. 
Jn Jan Jan 



2 Si 1 



zi; || | + / \Vw\ 2 dx + m / |tt)| 2 ds < c||u)|| 2 + / |/3 2 — /3i||tii||u;|<is. 
Jn Jan Jan 



It follows from /3j = min (max(— ^p-, m), M) , z = 1, 2, that 



!/?2 - /?i| < d u 2<£2 - ui<pi\ < ^(\u 2 ip\ + \wipx\). 
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Then we get 



- — |M| 2 + / \V ic\-<l:r + in j ir\-<ls 
1 



'n Jan 

|2 



2 



< c||«;|| 2 + — / \u2f2 — ui(f\\ \ui\ \w\ds 



< c \\ w \\2 + t; [ (\u 2 (f\ + \wipi\) \ux\ \w\ds 
Ian 



2 

< c|M|! + ^Wu^^^WoohWmd^Wwh^dn) + g Ho© llvi lloo||«7|li^(on) 
and, using Young's inequality, 

i^Hli + J Q \Vw\ 2 dx + fm- illmlloollyjilloo-cj HHIi^an) < c \\w\\l + c\W\\h(dny ( 22 ) 

On the other hand, using the adjoint system, we have 

dip 2 



dt 

and 

dip 
~ ~di 

where 



A<£ 2 = h 2 {ui,u 2 ,fi, <f2)(<P2 ~ <fi) (23) 



- - Atpx = hi(ux,U2,<pi,(p2)(<P2 - (pi) , (24) 



. f \ ( 2\J n f(u 1 )i Pl dx Xf(ui)(pi , i\ „ v 



and 



i, , \ ( 2A/ /(n 2 )y9 2 dx Af(n 2 )y? 2 \ 

Note that fai(iti, 1*2, ¥>i, <P7i), h2(u%,U2, (pi, <p%) € L°°(Q). Subtracting (j23|) from ([23]) . we get 

0(V?2 - 93l) 



dt 



+ A(v? 2 - </?i) = {hi(ui,u 2 ,(pi,(p2) ~ h 2 (ui,u 2 ,(pi,(p2)) (*P2 - (pi) 



Multiplying the above equation by (p = (p 2 — <pi, using hypotheses, and the L°° estimates of 
ui,U2,tpi,(f2,hi,h2, we get 

9 \\(f2 ~ (Pi\\l ~ I ! |V(v? 2 - ¥i)\ 2 dx - [ ^-{(f 2 - (pi){(P2 - (pi)ds\ < c\\<p2 - villi- 



2 9t {J n Jandu 

Then 

1 d 



2dt 



\~ I \^ip\ z dx+ I (Pi(pi - P 2 (f2) {(f2 ~ (pijds < c\\<f 11 o 
an 
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and it follows that 

ld_ 

2dt 



I - / \V<p\ 2 dx+ / f3 2 \<p\ 2 ds + / (f3 2 - p^ipds < c\\ip\\l 
jq Jan Jan 



We have = max (min(m, — ^^-),M). Therefore, 



\@2 ~ Pl\ < \U2V2 ~ U\(px\ < -{\ip2M + 



Using the fact that m < f3 2 , we have 

< [ \X7ip\ 2 dx + c\\<p\\l + - [ \u 2 <f2 ~ ui<pi\ \ipiip\ds 
Ju 1 Jan 

< I \V(p\ 2 dx + c\\ip\\l + - I {\<p 2 w\ + \ui<p\) \fif\ds 
Jo, 1 Jan 

f 1 1 

< J \Vlf\ 2 dx + C\\<p\\l + -\\f2\\oo\\fl\\oo\\w\\ L 2 {m) \\(f\\ L 2 {dn) + -||^l||oo||Ml||oo||y|li2 (M) . 

Using again Young's inequality, we get 

ia_ 

2dt 



+ yn - -||vi||oo||«i||oo - cj \\<p\\w(an) ^ J q l V( ^l dx + c IMl2 + c IHL=(an) (25) 

and, from Poincare's inequality and the fact that the operator trace from H l (£l) to the boundary 
space L 2 (dQ) is linear and compact, we have from ()22[) and ()25[) that 



2 dt 



Mb + imi2J 



) + (m- ^HuiHoolkilloo - cj \\w\\ 2 L2{m) 
+yn ~ ^IMUNHoc - cj \\(p\\ 2 L 2 {dn) < |Vv?| 2 (ix+c||?i;||^+c||^||^+c||u;||| 2(an) +c||^|! 



2 

L 2 (dQ)- 



Then, 



{W w \\l + IMIi) + (rn- ^||ui||oo|Mloo - cj \\w\\h(au) 



+ |m- ^IMUIMU - cj \\f\\ 2 L 2 im) < c|HI| + 
and for m sufficiently large one has 

^(Hli + IMII) <c(Hli + IMII)- (26) 

Gronwall's inequality leads to \\w\\ 2 + \\tp\\ 2 < 0. Then u\ = u 2 and (p 2 = </>i, which gives the 
uniqueness of solutions to the optimality system and therefore the uniqueness of the optimal 
control, since we have the existence of an optimal control and corresponding state and adjoint, 
which satisfy the optimality system. This completes the proof of Theorem 15.21 □ 
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Remark 2 : The uniqueness of the optimal control can be obtained from 

|/?2 ~Pl\< ^\U2^>2 ~ Uifi\ < ^{\u 2 f \ + \w<pi\), 

since ip = w = 0. 



6 Numerical Example 

We now give a numerical example for a particular problem. We use a finite element approach 
based on the Galerkin method to obtain approximate steady state solutions of the optimality 
system in the one-dimensional case. The formulation of the finite element method is based on a 
variational formulation of the continuous optimality system. The optimality system is discretized 
by finite differences. We then obtain the following one-dimensional nonlocal thermistor problem: 

du . Af(-u) 

8T- A " = a, /M*.)" 0<x<1 ' (>0 ' 

subject to the boundary and initial conditions 

On 

— = -(3u on aoxfo.r), 

dx 



u(x,0)=0, 0<x<l. 

We divide the interval £1 = [0, 1] into ./V equal finite elements = xq < x\ < . . . < xjy = 1. 
Let (xj, Xj+i) be a partition of and Xj+% — Xj = h = jj the step length. By S we denote a 
basis of the usual pyramid functions: 

\x + {l-j) on [ Xj -i,Xj], 
-\x + (l+j) on [xj,x j+ i], 
otherwise. 

First, we write the problem in weak or variational form. We multiply the parabolic equation by 
Vj (for j fixed), integrate over (0, 1), and apply Green's formula on the left-hand side, to obtain 



f du , r „ „ , f du , XL 
/ -rr v i dx + I vuvvj ax — — — Vj as = — -7- 
Jn 9t J n J dn du (J a 



x In f( u ) v J dx 

/(«) dxy 



Using the boundary condition we get 



/ !T v j dx + f VuVvj dx + ( (3uvj ds 



x Jnf( u ) v J dx 
(J Q f(u)dx)i- 



(27) 



We now turn our attention to the solution of system (|27p by discretization with respect to the 
time variable. We introduce a time step r and time levels t n = nr, where n is a nonnegative 
integer, and denote by u n the approximation of u(t n ) to be determined. We use the backward 
Euler-Galerkin method, which is defined by replacing the time derivative in (|27p by a backward 
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difference - + ~ u ■ So the approximations u n+1 admit a unique representation, 



TV 

n+1 n+l 



i=-l 



where a? +1 are unknown real coefficients to be determined. Thus, 



/■ n n + 1 — u n f f 

/ ujdx+ / Vu n+1 V^dx + / (3u n+1 Vj ds 

Jn T Jci Jan 



The scheme may be stated in terms of the functions vf. find the coefficients in u n+1 

Eil-i such that 



N N 

V < +1 / Vi Vj dx + rV] / V^V^ dx + r / /3n n+1 ^ (is 

__ 1 is] i= _ x is] Jan 



In matrix notation, this may be expressed as (.A + rB) a n+1 = g n = g(nr), where 

A = (flij) with element = / 



5 = (pij) with 6y = / SJviVvj dx , 



and a n+l is the vector of unknowns (a" +1 )^ = _ 1 . Since the matrix A and B are Gram matrices, 
in particular they are positive definite and invertible. Thus, the above system of ordinary differ- 
ential equations has obviously a unique solution. We solve the system (j28H for each time level. 
Estimating each term of (|28p separately, we have: 



N 



v^j dx 



i=-l 

/v 



iv x 

E < +1 / 

:=-i J o 

- a]+l / Vj-xvj dx + a] +1 / V J dx + c^ 1 / 



/ v'*dx+ Vjdx \ +a™^l / vjv j+ idx. 

JXi-i Jxj J J x j 
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Using the expression of Vj-i, Vj and Vj + ±, we obtain 



i=-l Jn 



Similarly, we have 



_ 1 

^ a * y n dx 



i=-l 

A? 



i=-l 

n+1 p ' dvj-idvj , , n+1 H +1 fd Vj \ 2 , n+i f X3+1 dvj dv j+1 

=-^i r d X+ °¥- r'dx-^si r +i dx 

- _Io< n+1 + -r/ i+1 - i« n+1 
" h a i- 1 + h a i h a i +1 - 



On the other hand, 



and 







9^={0,l} 



1 

J 2 



u n+ V ^ 4 (/3n" +1 (1)^(1) + /?« n+1 (0)^(0)) 



i(/3a^+S(l) + /K + S(0)) 

'i/3a£ +1 ifj = 0, 
- < if j = 1 . . . N - 2, 

if j = JV — 1. 



Furthermore, 
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Using the boundary conditions, we have 

a n _\ l = a «+ 1 + (/ l /3 + l) a ™+ 1 , 
a n _ x = a\ + (hj3 + 1) a£, 



a N 



n+l _ „n+l 



n ~ n 

From the initial condition we get a® = a % = ®- Setting 

h t 2h 2t 

and using together (|28p - (|33p . we then get the following system of N—l linear algebraic equations: 
for j = 0, 



a(l + + 6 + ag +1 + 2aa? +1 = ^ (5 + ^) a£ + ^aj + 
for j = l,...,N-2, 



(f(a$) + f(a%)r 



for j = N — 1, 



2h 



Similarly, 

^ n+1 = £ 



AT 

+1. 



i=-l 

where are unknown real coefficients to be determined. The discretization of the boundary 
conditions with respect to <p looks as follows: 

^f = ^ +1 + (/!/? + ik +1 , 

u "+i - i 

Mjv "l + /3/ i ^- 1 ' 

If we set 

/i r 2/i 2r 
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Figure 1. The evolution of temperature u 



Figure 2. The control 



then the remaining discrete equations, approximating the optimality system, are as follows: 
for j = 0, 



c{l + hp) +d + 



2At/3/'K) 



2 I + 2cnl 



n+l 



(/(«Sr) + /K)) 



6 3 (/K) + /K)) 3 



for j = 1, . . . N - 2, 



rn n+1 -4- Hn n+1 -4- rn n+i — tYi — — //" — — n n — — //" 



,71+1 



2/! 



1) 



for j = JV - 1, 



Finally, we have the discretization of [3 as follows: 



/ / u n+l^n+l 
(3 n+1 = min f max ( m, 



,M 



(34) 



The numerical experiments are in agreement with the results of (jSidi Ammi and Torresll2008l ): 
we obtain stable steady-state (see Figure [T]) . With an initial guess for the value of the control, 
the consecutive values of /3 converge to the lower bound when time is small and to the upper 
bound when t is big (see Figure [2]) . 
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